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Abstract. We consider a class of particle systems generalizing the /3-Ensembles from ran- 
dom matrix theory. In these new ensembles, particles experience repulsion of power /3 > 
when getting close, which is the same as in the /3-Ensembles. For distances larger than zero, 
the interaction is allowed to differ from those present for random eigenvalues. We show 
that the local bulk correlations of the /3-Ensembles, universal in random matrix theory, also 
appear in these new ensembles. This result extends the bulk universality classes of random 
matrix theory and may lead to a better understanding of the occurrences of random matrix 
bulk statistics in several observations which have no obvious connection to random matrices. 
The present work is a generalization of [GV12] where a similar result was proved for /3 — 2. 



Random matrix theory is well-known for universality phenomena which means that many 
essentially different matrix distributions lead in the limit of growing dimension to the same 
spectral statistics. 

In the past 15 years or so, much progress has been made in proving universality of local 
spectral distributions, especially correlations between neighboring eigenvalues in the bulk of 
the spectrum and of the largest eigenvalues. It is known that there is a parameter, usually 
denoted (3, which determines the universality class of the ensemble. To explain this in more 
detail, define for any /3 > and a continuous function Q : R — > R of sufficient growth at 
infinity, the invariant /3-Ensemble Pn,q,/3 on R^ which is given by 



(With a slight abuse of notation, we will not distinguish between a measure and its density.) 
For (3 = 1,2,4, Pn,q,/3 is the eigenvalue distribution of a probability ensemble on the space 
of (N x N) matrices with real symmetric (/3 = 1), complex Hermitian (/3 = 2) or quater- 
nionic self-dual (/3 = 4) entries, respectively. The matrix distributions are invariant under 
orthogonal, unitary or symplectic conjugations, respectively, explaining the name "invariant 
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ensembles". For arbitrary (3, only for quadratic Q, Pn,q,p is known to be an eigenvalue 
distribution. 

It has been shown (see |GV12] for references) that the local spectral statistics in the bulk or 
at the edges of the spectrum do in many cases not depend on Q or, in other terms, invariant 
ensembles with different potentials Q but the same (3 have the same local statistics. It is 
also known that different values of (3 lead to different limiting (local) distributions. This 
is not surprising as the interaction term Yl^j \ xi — Xjf has a strong effect on neighboring 

eigenvalues whereas e~ N '^i= 1 ®( Xj ^ just confines all eigenvalues independently into a compact 
interval. In the limit N — > oo, these two competing forces balance and produce a limiting 
measure of compact support. 

In |GV12| the question was addressed whether the interaction term Y\i<j I 

r could be 

changed without changing the local statistics. To this end, we introduced ensembles with 
density proportional to 



n^-'v •'•,)' v>: < ^-< : . (2) 

i<j 

where Q is a continuous function of sufficient growth at infinity compared to the continuous 
function ip : R — > [0, oo). The interaction potential tp fulfills 

(p(0) = 0, tp(t) > for t / and Jim ^ = c> for some (3 > 0, (3) 

or, in other terms, is the only zero of <p and it is of order (3. It has been conjectured in 
[GV12j that the bulk correlat ions for the ensembles ([2|) are the same as in the case ip(t) = \t\ , 
i.e. the same as for the invariant ensembles in random matrix theory. This was proved in 
GVI2J for (3 = 2 and a special class of functions <p and Q. In the present work, we prove a 
similar result for arbitrary /? > 0. This shows that the local bulk correlations (at least in the 
considered cases) merely depend on the repulsion exponent j3 and not on the interaction of 
eigenvalues at distances larger than 0. 

We believe that these results may lead to an explanation for the occurrence of random matrix 
bulk statistics in a number of seemingly unrelated observations in real world and science (see 
GV12J for references). Spacings between cars in different situations were found to be fitted 
well by the universal spacing statistics from random matrix theory ((3 = 1 for parking along 
one-way streets, (3 = 2 along two-way streets, (3 = 4 for waiting in front of traffic signals). 
Also spacings between perching birds and between bus arrival times at stops in certain cities 
seem to obey ((3 = 2) random matrix spacing statistics. Gaps between the Riemann zeta 
function on the critical line are another famous example from mathematics (also (3 = 2). In 
all these observations, a natural repulsion between consecutive quantities is present. 
Furthermore, the ensemble ([2]) does not seem to have a spectral interpretation which makes 
our findings a first step in proving universality of random matrix bulk distributions for more 
general particle systems. 

To state our main results, we first rewrite the ensemble ([2]). Let h be a continuous even 
function which is bounded below. Let Q be a continuous even function of sufficient growth 
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at infinity. By P^ g ^ we will denote the probability density on ~R N defined by 

1 N 
p n,qA x ) := ^h - xjf exp{-N^Q{ Xj ) -YsHxi ~ xj)}, (4) 

N,Q,/3 i<j j=\ i<j 

where Z% q p denotes the normalizing constant. The density -P/vq/3 can a l so be written in 
the form © with cp(t) := \tf exp{-/i(i)}. 

Furthermore, let for a probability density P/v on R^ and k = 1, 2, . . . , 

p N (ti, . . . ,t k ) := \ P N (ti,...,t k ,x k+ i,...,x N )dx k+ i...dx N 

jM. N ~ k 

denote the A;— th correlation function of P/v- The correlation functions are the marginal den- 
sities. The measure /?jv(ti, . . . ,t k ) dt\ . . . dt k is called A;— th correlation measure. Denote by 
p^ k q p the fc-th correlation function of P^ q ^ and by p k N q ^ the k-th. correlation function of 
Pn,q,i3 from (pTJ) . Universality of ensembles is usually defined by universality of their correla- 
tion functions or measures as many interesting statistics of the ensembles can be expressed 
in terms of correlation functions. Finally, introduce for a twice differentiable convex function 
Q the quantity a Q := inf tgR 

The following theorem deals with the global or macroscopic behavior of the ensemble P-fa q p. 

Theorem 1. Let h be a real analytic and even Schwartz function. Then there exists a 
constant a h > such that for all real analytic, strongly convex and even Q with oq > a h , 
the following holds: 

The first correlation measure Pj^q a converges weakly to a compactly supported probability 
measure plq a which has a non-zero and continuous density on the interior of its support. 
Weak convergence means that for any bounded and continuous f : R — > R ; we have 



lim / f(t)pf Qp (t)dt= / f(t)p h Q ^{t)dt. 



N- 

Remark 2. 

• In general, Pq a depends on h, i.e. changing the interaction term has an influence on 
the (limiting) global density of the particles. 

• If h is positive semi-definite, then a h in Theorem [T] may be explicitly chosen as a h = 
sup teR -/i"(t). 

• For k = 2, 3, . . . , the k— th correlation measure converges weakly to the /c-fold product 

p)® k '. This has been shown in [GV12] for (3 = 2 but the same proof goes through 
for arbitrary j3 > 0. As a proof of this statement would lengthen the presentation, we 
will not pursue it here. 

• Note that the dependence of pg^onf] can be eliminated if the prefactor (3 is put in 
front of Q and h. 

The next theorem states the local universality in the bulk. We use the notion of universality 
by Erdos, Yau et al. (see e.g. |EY12| and the references therein). Let G be the Gaussian 
potential G(t) := x 2 and recall that the corresponding limiting measure PQ,p is the semicircle 
distribution (with a certain variance depending on (3). Recall that under mild assumptions on 
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Q, there is a measure fiQ 3 of compact support which is the weak limit of the first correlation 
measure of Pn,q,/3- Consider the scaled correlation functions 



h,k 

,h f r ,\k P N,Q,(: 



a + 



, a + 



(5) 



where a is a point with fj,Q J a) > and t\, . . . ,t k are contained in an ^-independent com- 
pact interval. Under this scaling, the local density around a will be asymptotically one, in 
particular independent of a. For N — > 00, h = and Q = G, the limit of (J3j> exists and has 
been described in terms of a stochastic process in |VV 09j . As for general f3 no nice formula 
for this limit is known, we state the following theorem as universality result, comparing the 

with those of the Gaussian /3-Ensemble Pn,g,p- 



local correlations of Pjy 



Theorem 3. Let h and Q satisfy the assumptions of Theorem^ Let < £ < 1/2 and set 
sn '■= N~ 1+ £. Then for k = 1,2, ... , we have for any a in the interior of the support of 
Hq g, any a! in the interior of the support of the semicircle law fJ-c,/3 an d any smooth function 
f : M k — y M with compact support 



Jim I f(h,...,t k ) 



1 



a—spf 



h,k I 

■Pn,Q,P u + 



1 



a'+s N 

a'-s N VGA a ') 
= 0. 



KkPN,G,p u + 



tl 



N^ QB {a) 



,u + 



t-k 



tk 




du 




2s N 


du 
2s N 




. . dt k 



Remark 4. 

• If the inner integrations were not present, the convergence in Theorem [3] would be vague 
convergence of the scaled correlation measures. Here an additional small (uniform) 
average around the points a and a' is performed. 

• If h is positive semi-definite, then a h in Theorem [3] may be explicitly chosen as a h = 
supt eR -/i"(t). 

• The choice of the Gaussian /3-Ensemble Pn,g,/3 is just f° r definiteness, in fact any other 
ensemble belonging to the same universality class could be chosen. So far, these are 
known to be basically all Pn,q,/3 with the same f3 and real analytic Q which leads to a 
limiting measure jJ-n a of connected support [BEY12J. 



These results should be compared to those of }GV12j . There we could show for j3 = 2 under 
the same conditions on Q and h a much stronger type of convergence as in Theorem [3j We 
proved in [GV12j 



lim 



N^oo ^ h QJ3=2 {a) 



1 h,k I ti 

k PN,Q,8=2 a + Ml h 



^Q,/3= 2 («) 



det 



sin (7r(tj - tj)) 
n(U - tj) 



l<i,j<k 



, a + 



tk 
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uniformly in t±, . . . , tj~ from any compact subset of M and uniformly in the point a from any 
compact proper subset of the support of fig p_ 2 . This locally uniform convergence of the 
marginal densities was inherited from strong results on universality of unitary invariant (i.e. 
(3 = 2) ensembles (cf. [LL08] ) . In order to apply these results, we developed a method to 
express the correlation functions of the model Pjy q 0=2 as a probabilistic mixture of unitary 
invariant ensembles with potential V+f/N, where V was fixed and / was random. However, 
this representation was only possible for negative semi-definite h and an argument involving 
complex analysis had to be used to extend the universality for more general h. 
So far, the local relaxation flow approach due to Erdos, Schlein and Yau (refined by others) 
[ESYllj and applied to /3-Ensembles by Bourgade, Erdos and Yau [BEY111 IBEY12] is the 
only method for showing bulk universality for general /3-Ensembles. A remark on some crucial 
points of this method is included in Section 4. Their approach actually addresses universality 
of gap distributions which implies the weaker form of universality of the correlation measures 
as stated in Theorem El As we use their method, we obtain the same form of convergence. 
If other sufficiently general universality results on /3-Ensembles yielding stronger types of 
convergence were available, the method of [GV12j could be used to prove Theorem [3] with 
stronger forms of convergence. One advantage of the local relaxation flow approach is the 
possibility to compare local statistics of eigenvalue ensembles and other, not necessary spectral 
ensembles, directly. This allows us to give a short proof of Theorem [3l 

Similar in both |GV12] and the present work is the identification of the global behavior of 
Pjg o, e.g. Theorem [TJ The limiting measure is identified in Section 2 as a (unique) solution 
of a certain recursive equation. Using this measure, from Q and h a new potential V is defined 
and to P^Qp an invariant ensemble Pn,v,p is associated for which Pfa q a ~ exp{W}P/v,v,/3 
holds. In Section 3, a concentration of measure result is shown for U, furthermore Theorem [1] 
is proved. Section 4 contains the proof of Theorem [3] via the local relaxation flow approach. 
We finish with an appendix on some results about equilibrium measures which are needed in 
Section 2. 



2. The associated invariant ensemble 

In this section we determine the limiting measure for our particle system and use this to 
construct an ensemble of eigenvalues with the same global and local behaviour. This is 
analogous to [GV12] . In order to keep the presentation on the new features of Pjy q a self- 
contained, we include it here. 

Let f3 > 0, h be a continuous even function, Q a strictly convex symmetric function and 
assume that 

P^(x) := * JJ te-^fe-^Qto)-^***-**), ( 6 ) 

A N,Q,/3 l<i<j<N 
defines the density of a probability measure on H. , where 

Z h NO0 := [ TT \ Xi - Xj f e- N ^f=^'^ h ^-^dx 

' ' 7mjv i<ij<iv 



6 



denotes the normalizing constant. We will use the notation 

U(s):= f f(t-s)d(i(t), f^:= [ [ f(t-s)dn(t)dv(s) (7) 



for an even function / : M — > M of sufficient integrability. For the statement of the next 
lemma, M.\ will denote the set of compactly supported (Borel) probability measures on M. 
Recall that the unique minimizer of the functional 

IqAv) ■= f Q(t)Mt) + ^ f j \o g \s-t\- l d^{s)dii{t) 

is called equilibrium measure to the external field Q and (3 > 0. 

Lemma 5. Let h : R — > R be even, twice differentiable, bounded and such that h"{i) > —olq 
for allt. Define Th$ '■ M\ — ^ ^/i,/3(m) as the equilibrium measure to the external field 
t h-> Q(t) + hfj,(t) and ft. 

Then T^p has a fixed point, i.e. there exists a probability measure /j-q^ which is the equilib- 
rium measure to the external field 1 1— > Q(t) + J h(t — s)d//g pi s )- 

Proof. Without loss of generality we only consider the standard case (3 = 2 and omit the 
index /3, the general case can be reduced to this case by considering Q//3 and h/(3. 
Recall Schauder's Fixed Point Theorem which states that each continuous mapping T : C — > C 
of a compact, convex and non-empty subset C of a Hausdorff topological vector space has a 
fixed point. 

Consider the Hausdorff topological vector space Ai{K) of all signed finite Borel measures on 
some compact interval K of R, equipped with the topology of vague convergence. The subset 
Ail(K) of all symmetric Borel probability measures on K is non-empty, convex and compact 
by Helly's Selection Theorem. 

The first step in the proof is to show that it follows from h"(t) > — oq and the boundedness 
of h, that the support of the equilibrium measure to the external field Q(t) + h^(t) is for all 
[i included in a compact interval of R. By Theorem IA.6I (in the appendix), the support of 
the equilibrium measure for Q(t) + h^(t) is the smallest compact set K of positive capacity 
maximizing the functional 

K i y F Q+h ^(K) = logcap(K) - 2 / Q(t)du K (t) - 2 / h^{t)du K {t) (8) 
= F Q {K)-2 j h^{t)duj K {t). 

Choosing for K the support supp/i<g of the equilibrium measure [Iq, we get using \h^\ < ||/i||oo 
the simple inequality 

F Q+h„(.snppn Q ) > Fq(supp^q) - 2\\h\\ 00 G R. (9) 

It is easy to see that we have (h^)" = {h")^. By the condition h"(t) > — oq, Q(t) + h^{t) 
is convex for each compactly supported [i. By Theorem IA.6| the support of Th(fJ,) is a 
symmetric interval, say [—1^,1^]. Using Lemma IA.1[ we can rewrite ([8]) for an arbitrary 
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symmetric interval [—1,1] as 

FQ+h^H, I]) = Iog(i/2) - 2 J l _ Q(t)—^==dt - 2 £ h^t)—j^==dt. (10) 

Since Q is strictly convex and symmetric, we have Q(t) > aQt 2 + C for some CsK and (|10|) 
implies (using that the variance of is ^ 2 /2) the inequality 

*o+v(H. W < i°g(V2) - <*Q* 2 - C + 2||fc|| 00 , (11) 
which holds for any [i. From ([9]) and (lllh we see that 

^Q+fc^supp/iQ) > ^Q+^Q-M]), 

for all I > L for some /z-independent L. Hence L for all compactly supported 
Thus Tft maps the set M\(K) into itself, if X is chosen large enough. It remains to show 
continuity. We will show that Th maps converging sequences to converging sequences. Let 
C M}{K) be a sequence converging vaguely, or equivalently, weakly to a probability 
measure \x. Denote T/i(/z n ) ='■ v n . Then the sequence of external fields V n (t) := Q(t) + h^ n (t) 
converges pointwise to V(t) := Q(t) + h^(t). Since by Theorem I A. 41 the equilibrium measure 
does not depend on values of the external field outside of its support, we can assume this 
convergence to be uniform. Indeed, as h' is bounded on the compact set K by some constant 
C, we have \h'\ < C. It follows that the sequence {h^ n ) n is uniformly Lipschitz and 
hence equicontinuous. Thus the sequence (V n ) n is also equicontinuous. Since their domain 
is a compact and V n converges pointwise, the equicontinuity implies uniform convergence by 
Arzela-Ascoli's Theorem. 

Because all v n are supported on the same compact set, it follows that {y n ) n is tight and hence 
has a weakly converging subsequence [y nm ) m . We will prove that this limit measure, say u', 
is in fact v = Th(fi), the measure belonging the external field V, and does not depend on the 
particular subsequence. It follows that the sequence (v n ) n converges to v weakly. From the 
uniform convergence of V n towards V it follows by Theorem IA.51 1 . that 

TT^is) = J log\t - s]- 1 dv nm (t) 

converges uniformly (on C) towards U u {s) := J log \t — sp 1 dv(t). On the other hand, by 
Theorem I A. 51 we have for all s G C except a set of zero capacity 

lim U Vn ™(s) = U v '(s)= I \og\t- sr 1 dv'(t), 

m— >oo J 

which yields U u {s) = U v (s) almost everywhere on C and hence v = v' by Theorem IA.51 
implying that the sequence (u n ) n converges weakly to v. □ 

Remark 6 (Uniqueness). Uniqueness of the fixed point will follow for the class of ensem- 
bles from Theorem [TJ For those ensembles we will prove that the first correlation measure 
converges weakly to any fixed point, thereby showing uniqueness. 

We will now use Lemma [5] to construct an associated invariant ensemble. Let h be as in 
Lemma [5j Choose a fixed point fJq ^ as in Lemma [H We will stick to this measure from now 



on and write \i instead of fig p. Using notation (JJJ), we make the Hoeffding type decomposition 



N 2 N N 1 N 

j=l i,j=l 



T 2 »r N 



■—h^-—h(0) + N22ht,(xj)-U(x), where (12) 
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1 N 

U{x) := -- ( - Xj) - + ^(ac/) - VI)' ( 13 ) 

Now we can rewrite Pjy q s as 

^(s) = _J_ J] \ Xi - Xj fe- N ^ v ^*\ (14) 
Z W." i<i<J<iv 
where we defined the external field 

V(t) := Q{t) + h^t) 

and absorbed the constant exp{ — (N 2 /2)h fMfl — (N/2)h(0)} into the new normalizing constant 
Zn,v,/3,u- We will often use this representation of Pjy q «• The proofs of Theorems CD and [3] 
rely on comparison of Pjy n g with the invariant ensemble 

Pn,vA*)=t L - ]I \x t -x/e- N ^ v ^. (15) 

ZN > V >f> l<i<j<N 

3. Concentration of U 

This section is similar to Section 4 in [GV12j except for the proof of Theorem [TJ As several 
arguments of this section are needed later, we include it in the presentation. 
A key tool will be the following well-known concentration of measure inequality (AGZ10, 
Section 4.4]). 

Theorem 7. Let Q be two times differentiable with Q" > c > 0. Then we have for any 
Lipschitz function f and any e > 

N N 



cs2 

p N,Q^{\J2f(xj)-^N,Q,f3^2f(xj)\ >e) <2exp{- } and 

j=l j= i \ J \c 



- - e 2 \f\ 2 

E N>Q>g exp {e{ Y, ffrj) ~ ^N,Q,p £ /M } < ex P { 2c ^ 
j'=i i=i 
where we denote the Lipschitz constant of f by \f\jr- 

The following is a special case of a result in [Shell] (see also |KS10| ). By \xq^q we will denote 
the equilibrium measure to Q (and /3). 



Proposition 8. Let Q be a convex external field on R which is real analytic in a neighborhood 
of supp(/iQ )( g). Let f be a function whose third derivative is bounded on a neighborhood of 
supp(/xg i/3 ). Then 



N 

\E N , Q ,pJ2f( x j)- N / f d VQM <C(ll/lloo + ||/ (3) ||oo), 
i=i 3 

where C does not depend on N or f and \\ ■ ||oo denotes the bound on the neighborhood of 

SUPPOUQ,/}). 

Theorem [7] and Proposition [8] yield immediately 

Corollary 9. Let Q be a real analytic external field with Q" > c > 0. Then for any Lipschitz 
function f whose third derivative is bounded on a neighborhood of supp(//Q j( g) ; we have for 
any e > 

N r f 2 \f\ 2 

E^exp{e(^/(^)-A r / f(t)dfx Q ,^t)) } < exp {-^ + eC(\\f\\oo + ll/ (3) IU)} • 
i=i 3 ° 

The idea is to reduce concentration of the bivariate statistic IA to concentration of linear 
statistics. To this end, we give an alternative representation of IA using Fourier techniques. 
A similar idea is used in [LP08j . 

Lemma 10. We have 



U{x) = -= I iiN(t,x) h(t)dt, where 



2V2vr 

N f N ~ 1 f 

U]\[(t,x) := y cos(txj) — N cos(ts)dfi(s) + V—l } sin(txj), h(t) := / e 

j^l J ~y v27r Jr 

Proof. Recall from (|13p that 

1 N 

U{x) = --( ^ - Xj) - [/i M (xj) + /^(xj) - ). Note that 

^ frfo - x k ) = -±= [ J2e i(Xj - Xk)t Ht)dt = -±= [ \u N (t,x)\ 2 h(t)dt, 
2 j,k 2V2lT 3 j,k 2V27T 3 

with UN{t,x) := J2f=i eltXj - Writing uj<f(t, x) := U7v(t,x) — N f e lts d/j,(s), it is not hard to 
check that 

1 f 2 ^ 

U(x) = =/ u N (t,x) h(t)dt. (16) 



h(s)ds. 



2V27T 

□ 

A trivial but useful observation is 

EjV )Q)/3 /(z) = (Z N y^/Z Ny ^ u )W, N y^f(x)e U ^. 
The next proposition establishes concentration of IA. 
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Proposition 11. If the constant uq is large enough, then there exist constants C\, C2 > 
such that for all N 

< Ci < Z N y i p tU /Z N y }j3 = Ejv^exp {U(x)} < C 2 . 
Proof. By Jensen's inequality we get 

Ejv,v,/9 exp {U(x)} > exp 
Using Lemma[10]we show that the expectation of hi is bounded in N . Fubini's Theorem gives 

If o 2 - 

— E NtVj pU(x) = — = / E N yp UN(t,x) h(t)dt 
2v27r J 

1 



2V2tt 



N „ N 

EjsiyA 



cos(txj) — N / cos(ts)e?//(s)| 2 + Ejv.y^l^^ sin(txj)| 2 J h(t)dt. 
i=i J 3=1 ) 



By Corollary [H the terms in the parenthesis are bounded by a polynomial in t, as 
|cos(t-)|£ , |sin(t-)|£ < t and ||cos(t-)*- 3 ' 1 ||oo, ||sin(t-)^ 3 ^ ||oo < Ct 3 . Hence, h being a Schwartz 
function, we have E^y^U{x) > — C for some C > 0. Thus the lower bound follows choosing 
d :=exp(-C")- 

For the upper bound recall that since h is even, h is real- valued. Define h+{y) := max{0, h(y)} 
and h-(y) := max{0, —h(y)} such that h = h + — h—. For h- = 0, which corresponds to the 
case of a positive semi-definite h, there is nothing to prove, so assume that /i_ 7^ 0. 

^ 1 /2 

Introducing iT_ := f/i-) > 0, Jensen's inequality and Tonelli's Theorem give 



E N y 



;/3 exp| - (2v / 2^) _1 y h(i)|«jv(i,x)| 2 di} < E N>VtP exp j^v 7 ^) -1 ^ #-(t) 2 |«Jv(*, x)| 2 dt} 

exp{(2v / 2^)- 1 ||i/-|| Li y (!T_(t)/||ff_|| Ll )fl-_(t)|u jV (t,a;)| 2 (it} 

< J (^_(0/||^-|| L1 )E^exp{(2v / 27)- 1 ||^-|| Ll F-(t)|^(t,x)| 2 }df. (17) 
Abbreviating := (2v / 27r) _1 ||i?_|| i . 1 and using the Cauchy-Schwarz inequality, we find 

E N yp exp {K h H- (t) \u N (t,x)\ 2 } (18) 

< fij^exp ^2K h H-{t) J^cob^) - iV / cos(ts)dfj,(s) } (19) 

iV 2 

x exp { 2i^#_ (i) I ^ sin(t Xj ) | } . (20) 

i=i 

Considering only (|19p . we have by Corollary [9] for any e > 0, 

N 

E7v,y,/3exp je • y / 2K h H_(t) ( ^ cos(ixj) - iV / cos(ts)dfi(s)) } 
< exp {e 2 • 2K h H_(t)t 2 (2a v )- 1 + e^/2K h H_(t)C(l + i 3 )}, (21) 
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where ay '■= mint V"(t) > an d C does not depend on t or N. For oq large enough (hence 
ay large enough), we have 2K h H_(t)t 2 (2ay) -1 < 1/4 for all t. Since i3_(i) = h l l 2 {t) is 
decaying rapidly, y '2KhH ~{t)C '(1 + i 3 ) is bounded in i. Summarizing, if ag is large enough, 
we can bound ([2T]) by 

exp{ce 2 + eC} 

with < c < 1/4 and c, C do not depend on iV or t. Thus we have 
^jv,v,/3 exp ^2-K"faJ3"_(t) y~^cos(tej) — N I cos(ts)dfj,(s) j 



3=1 

N 



— / exp < e ■ ^2K h H_(t)( V cos(txj) — N cos(ts)dfj,(s)) \ exp{-e 2 /4}cfe 
2tt I ~— J ' J ) 

< — exp{ce 2 +eC}exp{-e 2 /4}de < C for some C . 

We conclude that (|18p is bounded in TV as long as ag is large enough. Finally, it follows from 
(fP7j) that 

K7v,y,/3 exp {— J h(t)\u^(t, x)\ 2 dt} < C 

for some constant C > independent of N. This proves the upper bound. □ 

Remark 12. The proof of Proposition [11] actually shows that for each A > there is a 
threshold a h (\) > and constants C\,Ci (depending on A and a h ) such that 

< d < E Ny ^exp{XU(x)} < C 2 , if a Q > a h (X). 

Proof of Theorem [7J It remains to show that the fixed point fx, whose existence was obtained 
in Lemma [5l is unique and indeed the limit of the first correlation function. We consider a 
Lipschitz function / : M — > M. with three continuous derivatives and estimate for any e > 

N 



PN, Q d\ N ~ l Y.tt X i) ~ J fM > £ ) = (ZN,V,/3/Z N y^ U )^N,V,^ X h {lN _ lj: N =i 



Lif(xj)-ffdl*\>eY 

3 = 1 

By Holder's inequality and Remark 1 12} we have 

N N 

^Wl^E-^)- fdv\>e)<C{P Ny A\N- 1 Y,f(*3)- I fd^\>e)) c 

3=1 3=1 

for some c,C > 0. By Corollary [SJ this last probability converges for any e > to 
exponentially fast as N — > oo. We conclude with Lebesgue's Theorem that 

lim E^liV- 1 ^/^)- / fdn\=0 and hence lim E^JNT 1 £ f( Xj ) = / fdfi. 

3=1 J 3=1 J 

As convergence for smooth Lipschitz functions determines weak convergence, the weak con- 
vergence of the first correlation measure follows. As the limit of weak convergence is unique, 
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this shows uniqueness of the fixed point in Lemma [5j The existence and positivity of the 
continuous density of is clear by Theorem IA.7I as V is real-analytic and strictly convex. □ 

4. Proof of Theorem [3] 

In this section we use the local relaxation flow approach developed by Erdos, Yau, Schlein et. 
al. to establish universality of the local bulk correlations. First we introduce some notation 
from [BEYllj . 

Let k be fixed. Let G : M. k — > R be a smooth function with compact support and m = 
(mi, . . . , mfc) with rrij being positive integers. Define 

Gi^ m (-c) • — G{N {xi Xi-t- mi ),..., A^(xj_|_ mfe _ 1 Xj-|_ mfe )) • 

The Dirichlet form of a smooth test function / : M. N — > R w.r.t. a probability measure du 
on R N is defined as 

N 



DM) : =^E/ {d Xj ffdu. 



3=1 

Let / be a probability density function w.r.t. du. The (relative) entropy of / w.r.t. du is 
defined as 

S U (J) ■= J flogfdu. 
We will use the following general theorem. 

Proposition 13. [B EYlll Lemma 5.9] LetG : M. k — > R be bounded and of compact support. 
Let du be a probability measure on {x : x\ < X2 < • • • < xn} C R n given by 

du = ie-^^W, H(x) = H (x) - ^ E log \x t - x { \ 

i<j 

with the property that V 2 %o > t _1 holds for some positive constant r. Let qdu be another 
probability measure with smooth density q. Let J C {1, 2, . . . , N — — 1} be a set of indices. 
Then for any E\ > we have 



-c7V e i 



In our application we will choose du = Pn,v,/3 and g = {Zj^yp/Zj^y^^u) exp{U}. If oq is 
large enough, then V is strongly convex, hence r = 1/ay. By the symmetry of Pjy-g p and 
Pn,v,p-, it is equivalent to restrict the measure to the simplex {x : x\ < xi < • • • < xjy} and 
multiply by iV! . From [EY121 Theorem 2.3] we have that 

SM) < cduVq)- 

It is thus sufficient to prove that D ul {^fq) is bounded in N as J will be chosen such that 
I J I ~ N in order to identify the bulk correlations. 



13 



Remark 14 (On the local relaxation flow approach). To briefly explain the essence of this 
method due to Erdos, Schlein, Yau and others (see e.g. |EY12] for references and a complete 
review), let us consider two measures as in Proposition 1131 duo and qduo and their statistics 
J gduo and J gqduo for some test function g. Assume that one can define a Markov process 
on M. N in terms of the Dirichlet form D w (or the formal generator Ln := t^t A — |(V'H)V), 
having duo as stationary distribution. Assume that the process has the initial distribution 
qduj and denote the evolution of the density w.r.t. dw by (ft)t>o, fo = Q, foo = 1- Then one 
can write 

/ 9 qdu- J 9 duo=(J gqduo- J gf t duo) + (J gf t duo- J gduo), 

which corresponds to running the process up to time t. If the process is ergodic and the time 
t is large enough, j g ftdu will be close to the equilibrium J gdoj. If this t is still "small", 
i.e. the convergence to the stationary distribution is fast, then the distance between J gqdu 
and J g ftdu should be not too big. These distances are measured in terms of Dirichlet 
form and entropy of duo. These estimates are due to the Bakry-Emery method which yields 
ergodicity or relaxation making use of the strict convexity of the Hamiltonian, i.e. of the 
bound V 2 %o > t . It turns out that the constant r is the time scale for the relaxation 
to equilibrium, meaning that e.g. Su(ft) < e~*/ T 5' w (/o). Here we tacitly used that the 
logarithmic part of the Hamiltonian % is convex, therefore does not increase the relaxation 
time. However, one crucial observation is that from the trivial bound 

(v, V 2 H(x)v) > -\\vf + 4 E ^"^2 

one can infer that the relaxation is much faster in the directions (uj — Vj) provided that Xi 
and Xj are close. Indeed, the mean distance between neighboring eigenvalues is of order 1/JV, 
hence the convexity bound for the Hamiltonian should be locally of order N, therefore yielding 
a time to the local equilibrium of order 1 /N whereas the time to the global equilibrium is of 
order 1. This informal reasoning can be captured by choosing test functions like Gi )Jn which 
depend only on eigenvalue differences in the local scaling (i.e. multiplied by N) and vanish 
whenever two eigenvalues are not close to each other. By exploiting these features of Gi iTn 
and some estimates, one arrives at Proposition [T3j For arbitrary test functions <?, one would 
get basically the same estimate except for the quantity | J\ ~ N which divides D u (y/q). 
One problem with this idea is that the existence of the process associated to the Dirichlet 
form is not clear for /3 G (0, 1). For /3 > 1, the repulsion is strong enough to prevent 
collision between the eigenvalues but for f3 < 1 the probability of explosion is positive. This 
problem was overcome in [EKYY12J by smoothing the singular logarithmic term and using 
the approach above for the corresponding process. 

In most applications the method sketched above was not directly used to prove bulk univer- 
sality because the quantity D L0 {^fq)/\J\ could not be estimated to decay as — > oo. To 
overcome this difficulty, duo was modified by adding another strongly convex potential to the 
Hamiltonian which confined the eigenvalues Xj near their classical positions 7j which are de- 
fined by the relation j = N d/j,, where \i is the equilibrium measure of the ensemble. This 
additional convexification was shown to slow down the relaxation time to the global equilib- 
rium t such that t ~ N~ £ which was then enough to prove bulk universality. To show that 
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the local statistics are not influenced by this modification of the ensemble, a strong bound on 
the rigidity of the eigenvalues is needed, which means informally speaking concentration of 
the eigenvalues around their classical positions. This strategy of proving bulk universality has 
proven very useful whenever strong bounds on the rigidity of the eigenvalues are available. 
Much effort in the works of Erdos, Yau et al. has been made to establish these bounds. 
In our application, no additional convexification of the Hamiltonian is needed as we can 
effectively estimate the Dirichlet form in our case. This is due to the fact that duj and qdu 
have the same global limit and we have concentration of IA under du = Pn,v,p- 

Proposition 15. Let D^y,p denote the Dirichlet form w.r.t. Pnvb an d 

q = {Z^y^/ZNy^pi) exp{£Y}. Then there is a constant C such that we have for oq large 

enough 

D N yAVQ)<C for all N. 

Proof. The ratio Zn,v,b/Zn,v,P,U ls bounded by Proposition [TT] and therefore negligible. We 
have by Holder's inequality for e > 

N N 

D N ,vAVQ) E ^ ( 9 *i exp{-^(x)}) 2 = C— ®N,V,8 exp{U(x)} {d Xf U(x)) 2 

1=1 i=i 

1 N 

< C(E N y P eM(l + e)U(x)}) 1/£ — £ {^ N y P \d Xl U{x) \*+Wy/<F+V. 



8N 
1=1 



Again by Proposition 1111 (Ejv.v^s exp{(l + e)U(x)}^ 1 ^ is bounded in N. In order to bound 
the second term, recall that 

i r N r N 

U(x) = — I (I y^cos{txj) - N I cos(ts)d/i(s)\ 2 + \y2sm(txj)\ 2 )h(t)dt. 

2V2ttJ j=l ' J j=1 



In the following we only treat the cosine term, the term involving the sine can be estimated 
analogously. We have 

N „ 



|2(e+l)/e 



V J cos(ts)dfi(s)\ 2{£+1)/£ \t\ 2{£+1)/£ \h{t)\dt (22) 

where the last inequality is derived by first applying the triangle inequality and then using 
Jensen's inequality. 

Taking now expectations, we get from Corollary [9] and the strong decay of h that the expec- 
tation of (|22p is bounded in N. This gives the claimed bound. □ 





N 




(|J^coe(ta?j) 












(^cos(tej) 








. N 


<_cj 


| cositxj) 
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Proof of Theorem [2 From Propositions[[3]and[l5]we have that the statistics fjr Yliej q p^i 
and tjt Y^ieJ^NjVtpGi,™ coincide in the limit N — > oo, as long as limjv->-oo = for some 
£\ > 0. It is a standard argument ( |ESYYf2| Section 7]) to infer from this that also the cor- 
relation measures of Pj^ q g and Pn,v,p coincide in the sense of Theorem [3j The universality 
of these correlation measures is the universality result [BEYlll Corollary 2.2] which precisely 
states that the correlation measures of -P/v,Qi,^ and Pn,q 2 ,P have the same limit (in the sense 
of Theorem [3]) for any real analytic and strongly convex Q\, Q2 with oq 1 , cxq 2 > 0. □ 

Appendix A. Equilibrium Measures with External Fields 

In this appendix, we recall some results about equilibrium measures, mainly from the book 
by Saff and Totik [ST97, Section 1.1]. The following can be found in [ST97, Section 1.1]. Let 
7W 1 (S) denote the set of Borel probability measures on a set E. Define for E C C compact 
the logarithmic energy of \i E Ai 1 (Ti) as 

I(jj):= I log \z - 1]' 1 dn(z)dn{t) (23) 



and the energy V of S by V := inf^g^i^) It turns out that V is finite or 00 and in 

the finite case there is a unique measure which minimizes ()23p . This measure is called 
equilibrium measure of S and the quantity cap(S) := e~ v is called capacity of E. For an 
arbitrary Borel set S we define the capacity of E as 

cap(E) := sup{cap(i^) : K C S compact}. 

Lemma A.l. If E = [—1,1], I > 0, then cap(E) = 1/2 and the equilibrium measure is the 
arcsine distribution with support [—1,1]: 

du^t) = —j^==dt, te[-l,l}. 

LdY, has mean and variance l 2 /2. 

Proof. See [5T971 Section LI]. □ 

Definition A. 2. Let Sclbe closed. Let Q : E — > [0, 00] satisfy 

a) Q is lower semicontinuous, 

b) Eo := {t £ E : Q(t) < 00} has positive capacity, 

c) if E is unbounded, then lim Q(t) — log \t\ = 00. 

|t|^oo,tes 

If Q satisfies these properties, we call it external field on E and W = e~® its corresponding 
weight function. 

Furthermore, define for fi £ A / 1 1 (E) the energy functional 

Iq(p)'= f Q(t)dn(t) + [ /log|s-tr 1 ^(s)(i/i(t). (24) 



Remark A. 3. In [ST97] the authors define the energy functional to be (in our notation) I2Q 
instead of Iq. It is more convenient for our purposes to use this definition. We note that 
under this change qualitative results from |ST97j remain the same but quantitative results 
involving Q have to be changed by a factor 2 or 1/2, respectively. 
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Iq(h) might be oo, but the following theorem holds. The support of a measure /i will be 
denoted as supp^). 

Theorem A. 4. Let Q be an external field on X. 

a) There is a unique probability measure /iq G A / 1 1 (S) with 

Iq(jjlq)= inf Jq(m). (25) 

b) yUQ /ias a compact support. 

c) Let Q be an external field on £ smc/i £/tcrf Q = Q on a compact set K with supp(/xg) C K 
and Q(t) = oo /or t ^ K. Then /ig = /ig. 

Proof. Statements 1) and 2) can be found in [ST971 Theorem 1.1.3], 3) follows from [ST97, 
Theorem 1.3.3] (also see remark on page 48 in [ST97J). □ 

/ig is called the equilibrium measure for Q. The next theorem summarizes properties of the 
logarithmic potential 



U^{z) := J \og\z-t\~ 1 d{i{t). 



Theorem A. 5. 

a) Let Q and Q be external fields on £ such that \Q — Q\ < e on S. JTten /or all z £ C 

\U^{z) - U^Q{z)\ < 2e. 

b) Let ii - C i fc compact and (fi> n )n be a sequence in M l {K) converging weakly to a 
probability measure [x. Then for a.e. zeC (w.r.t. the Lebesgue measure on C) 

liminf U^ n (z) = U^(z). 

n— >oo 

c) If fi and v are two compactly supported probability measures and their logarithmic 
potentials and U u coincide almost everywhere on C, then [i = v. 

Proof. Statement 1. is contained in ST97, Corollary 1.4.2], statement 2. is |ST97|, Theorem 
1.6.9] and assertion 3. [ST971 Corollary II.2.2]. □ 

We also need a characterization of the support of the equilibrium measure. 

Theorem A. 6. Let Q be an external field on S. 

a) For a compact set K of positive capacity define the functional 

F Q {K) := logcap(^) - 2 / Qdw K . 

For any compact K of positive capacity we have Fq(K) < Fq (supp(/iQ)). Furthermore, 
if K is compact and of positive capacity and such that Fq(K) = -Fq(supp(/zq)), then 
suppOq) C K. 

b) If Q is convex, then supp( / ug) is an interval. 

c) If Q is even, then supp(^g) is even. 

Proof. For statement 1. see [ST971 Theorem IV. 1.5], for statements 2. and 3. |ST97|, Theorem 
IV.1.10]. □ 
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Theorem A. 7. 

a) Let Q be an external field onS. If Q is finite on supp(^Q) and locally of class C l+£ for 
some e > (which means that Q is continuously differ entiable and the derivative Q' is 
Holder continuous with parameter e ), then \iq a continuous density on the interior of 
suppOg). 

b) If Q has two Lipschitz derivatives and is strictly convex, then supp(/ig) =: [a, b] and 
the density of [iq can be represented as 

dji(t) 



dt 



r(f)V(t-o)(6-t)l M] (i), 



(26) 



where r can be extended into an analytic function on a domain containing [a, b] and 
r(t) > fort 6 [a, b]. In particular, density is positive on (a,b). 



Proof. Statement 1. is |ST97} Theorem IV. 2. 5], for assertion 2. see e.g. the appendix of the 
paper by McLaughlin and Miller [MM08j . □ 
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